A posteriori error analysis of a quadratic finite volume method for nonlinear elliptic problems

In this article, we construct and analyze a residual-based a posteriori error estimator for a quadratic finite volume method (FVM) for solving nonlinear elliptic partial differential equations with homogeneous Dirichlet boundary conditions. We shall prove that the a posteriori error estimator yields the global upper and local lower bounds for the norm error of the FVM. So that the a posteriori error estimator is equivalent to the true error in a certain sense. Numerical experiments are performed to illustrate the theoretical results.     

Efficient Collocation Schemes for Singular Boundary Value Problems

We discuss an error estimation procedure for the global error of collocation schemes applied to solve singular boundary value problems with a singularity of the first kind. This a posteriori estimate of the global error was proposed by Stetter in 1978 and is based on the idea of Defect Correction, originally due to Zadunaisky. Here, we present a new, carefully designed modification of this error estimate which not only results in less computational work but also appears to perform satisfactorily for singular problems. We give a full analytical justification for the asymptotical correctness of the error estimate when it is applied to a general nonlinear regular problem. For the singular case, we are presently only able to provide computational evidence for the full convergence order, the related analysis is still work in progress. This global estimate is the basis for a grid selection routine in which the grid is modified with the aim to equidistribute the global error. This procedure yields meshes suitable for an efficient numerical solution. Most importantly, we observe that the grid is refined in a way reflecting only the behavior of the solution and remains unaffected by the unsmooth direction field close to the singular point.     

A likelihood-based method for haplotype association studies of case-control data with genotyping uncertainty

This paper discusses the associations between traits and haplotypes based on Fl (fluorescent intensity) data sets. We consider a clustering algorithm based on mixtures of t distributions to obtain all possible genotypes of each individual (i.e. "GenoSpec-trum"). We then propose a likelihood-based approach that incorporates the genotyping uncertainty to assessing the associations between traits and haplotypes through a haplotype-based logistic regression model. Simulation studies show that our likelihood-based method can reduce the impact induced by genotyping errors.     

棱镜自重变形对波面影响的研究

采用有限元方法分析了单棱镜在重力作用下三个典型位置的表面变形，对棱镜表面Φ425mm通光孔径轴向变形量进行Zernike多项式拟合，画出了相应的波差图，并确定设计方案.进一步选择Φ360 mm口径光束对双棱镜在几种典型相对位置下进行光线采样、追迹，计算得到不同情况下由于双棱镜自重变形引起的总波差.     

Effective magnetic properties of a composite material with circular conductive elements

Effective magnetic properties of a composite meta-material consisting of periodically arranged circular conductive elements are studied theoretically. A general expression for the effective bulk permeability is obtained with mutual effects and lattice ordering being taken into account. The resonance frequency of the permeability is found to be strongly dependent on the size and shape of the unit cell. Frequency dispersion of the permeability is studied with special attention paid to the frequency range, where negative values of the permeability are possible. Corresponding recommendations for optimisation of the meta-materials with negative permeability are made. The results are confirmed by numerical simulations of the finite structure behaviour in an external magnetic field. Received 19 April 2002 Published online 31 July 2002     

An efficient numerical scheme for precise time integration of a diffusion-dissolution/precipitation chemical system

A numerical scheme based on an operator splitting method and a dense output event location algorithm is proposed to integrate a diffusion-dissolution/precipitation chemical initial-boundary value problem with jumping nonlinearities. The numerical analysis of the scheme is carried out and it is proved to be of order 2 in time. This global order estimate is illustrated numerically on a test case.

     

关于AOR迭代法的研究

本文论证了严格对角占优矩阵之AOR法的误差估计式中的误差估计常数hγ,ω（0≤γ≤ω0）的最小值是h1,1.     

A dynamical test of phase transition order: New things in old places or old wine in new bottles

We first discuss nonlinear aspects of phase transition theory applied to a particular liquid crystal phase transition. A simple derivation is given to show how two coupled Goldstone modes (one appearing as gauge fluctuations of the ordered phase) can force a phase transition, against all expectations, to take place discontinuously (theory of Halperin, Lubensky, and Ma)-but the discontinuity may be immeasurably small. Then, we describe a new dynamical test of phase transition order, developed by Cladiset al., that turns out to be more sensitive than x-ray diffraction and adiabatic calorimetry. Quantitative data found by this new method are in excellent agreement with the measurements of adiabatic calorimetry and x-ray diffraction as well as expectations implicit in the predictions of HLM.This is the text of an after-banquet talk given at the CNLS Workshop on the Dynamics of Concentrated Systems.     

校正空间分块镜复位误差所需的最少自由度的确定

可收展式分块镜共相位成像技术是实现光学系统超大口径、超高分辨率成像的关键技术。分块镜复位误差会引起各子镜之间的共相位误差和整个拼接镜面的面形误差。合理确定校正分块镜复位误差所需的自由度具有非常重要的意义。计算机仿真分析了分块镜的复位误差对拼接镜面面形误差的影响。结果表明,所需的校正自由度的个数不仅与展开机构的复位精度有关,而且与分块镜镜面的几何参数有关。     

Discretization of the solutions to Poisson’s equation

Sharp estimates (in the power scale) are obtained for the discretization error in the solutions to Poisson’s equation whose right-hand side belongs to a Korobov class. Compared to the well-known Korobov estimate, the order is almost doubled and has an ultimate value in the power scale.